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ABSTRACT 

Frequency shifts, radiative decay rates, the Ohmic loss contribution to the nonradiative 
decay rates, fluorescence yield, and photobleaching of a two-level atom radiating any- 
where inside or outside a complex spherical nanoshell, i.e. a stratified sphere consisting 
of alternating silica and gold concentric spherical shells, are studied. The changes in the 
spectroscopic properties of an atom interacting with complex nanoshells are significantly 
enhanced, often more than two orders of magnitude, compared to the same atom in- 
teracting with a homogeneous dielectric sphere. The detected fluorescence intensity can 
be enhanced by 5 or more orders of magnitude. The changes strongly depend on the 
nanoshell parameters and the atom position. When an atom approaches a metal shell, 
decay rates are strongly enhanced yet fluorescence yield exhibits a well-known quenching. 
Rather contra-intuitively, the Ohmic loss contribution to the nonradiative decay rates for 
an atomic dipole within the silica core of larger nanoshells may be decreasing when the 
silica core - inner gold shell interface is approached. The quasistatic result that the radial 
frequency shift in a close proximity of a spherical shell interface is approximately twice 
as large as the tangential frequency shift appears to apply also for complex nanoshells. 
Significantly modified spectroscopic properties (see computer program freely available at 
|http: / / www.wav e-scattering.com ) can be observed in a broad band comprising all (non- 
resonant) optical and near-infrared wavelengths. 

PACS numbers: 78.67.Bf, 33.70.Jg, 32.70.Jz, 33.50.-j, 87.64.Ni, 87.64.Xx 
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1 Introduction 



Spectroscopic properties of an isolated atom, such as radiative and nonradiative decay 
rates, frequency shifts, and fluorescence yields are not inherent to the atom but charac- 
teristic of an atom coupled to a physical system. Indeed, it has been known for a long 
time that the presence of a small structure, cavity, or an interface can significantly im- 
pact the characteristic behaviour of a radiating system, irrespective if the emission takes 
place inside or adjacent to a material body [H 121 El- The physical origin of the decay of 
an excited atom state is the coupling of the atom to the vacuum electromagnetic field. 
A nearby presence of a material body modifies the vacuum electromagnetic field at the 
atom position. Consequently, the atom interacts with the modified vacuum electromag- 
netic field and it will exhibit different spectroscopic properties than the same atom in 
the absence of the material body. There is growing interest in the application of various 
systems that can significantly affect the vacuum electromagnetic modes. Such systems are 
currently of great interest in the fields of photonics and quantum electrodynamics. They 
have found widespread application in microcavity lasers, electroluminescent devices, and 
proposed photonic band-gap devices. There is also growing interest in the application 
of various aspects of the moleculesurface interaction in the field of medical diagnostics, 
particularly in the immuno-assay area, in which fluorescence-based techniques are widely 
used 0. Alternatively, for the applications in near- field optical microscopy, one is inter- 
ested in changes in the atom fluorescence properties induced by the presence of a nearby 
microscope tip. 

For the purpose of this article, the atom would mean any localized fluorescent dipole 
source, e.g., fluorescent organic group, rare earth atom, etc. The atom would be considered 
as a two-level system in the regime of weak coupling, within the domain of applicability 
of the linear response theory El IE] In the latter case, the quantum- mechanical 
description |Hj yields identical results to the classical description [HI CHI UH E21 
E31 Ej. The spectroscopic properties of the atom will be studied as a function of the 
atom position inside and outside stratified spheres. The theory of fluorescence properties 
of an atomic dipole has has been mostly investigated only in the case of a homogeneous 
sphere jH |H1 Ell HH EH EH EH] • Chew et al [TB] provided a formal solution to the problem 
of a dipole radiating in the presence of a multilayered sphere. However, their solution 
for the sphere with N concentric shells (the sphere core counts as shell number one) is 
written in terms of a 2N x 2N matrix and appears awkward and impractical for numerical 
calculations. Indeed, neither Chew nor anybody else have appeared to implement the 
Chew et al jTHj solution numerically. The main obstacles are that as N increases so do 
computer memory requirements to store the matrix and the time to carry out the matrix 
calculations, which increases as N 3 . Whenever a radiating dipole has been discussed 
interacting with a multicoated sphere, either the problem has been treated in a quasi- 
static approximation [3;, or the dipole has only been allow to radiate at the sphere origin 
pj EH EH EH EHj • The dipole position in the center of such a complex sphere considerably 
simplifies calculation as the inherently vector problem reduces to the scalar one involving 
scalar fields r ■ B(r) and r • D(r) ^71 EB]- Additionally, in most cases only the simplest 
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core-shell particles have been dealt with [7J HH1 1201 • Although Li et al [21] has provided 
a recursive formula for Green's function for the case of a multicoated sphere with an 
arbitrary number of concentric shells, an efficient numerical use of Li's formula, even 
for coinciding spatial arguments, requires to perform traces over the magnetic angular 
momentum number. Only very recently the limitation on the dipole position at the center 
of a general multilayered sphere has been removed and the traces over the magnetic angular 
momentum number in the scattering Green's function at coinciding spatial arguments, 
G(r, r, u), have been explicitly performed [12]. A complete description of the classical 
electromagnetic fields of a radiating electric dipole has been achieved outside and inside a 
multi-structured spherical particle. Electromagnetic fields have been determined anywhere 
in the space, and the time-averaged angular distribution of the radiated power, the time- 
averaged total radiated power, radiative and nonradiative (due to Ohmic losses) decay 
rates, frequency shifts have been calculated. Our recursive solution only employs 2x2 
transfer matrices and their ordered products and provides a fast and reliable algorithm 
which can easily be implemented numerically [T2"] . 

In the present article, the theory developed in Ref. ^2] will be applied to "nano- 
matryoshka" structures of Prodan et al [22], i-e., multilayered spheres consisting of alter- 
nating silica and gold concentric spherical shells. Such complex spheres have also been 
known as nanoshells. Current experimental colloidal techniques allow one to design a 
variety of multi-structured beads having a plurality of concentric shells with the core ra- 
dius from cca 1 nm till 1 fim and controlled shell thicknesses. For instance, metal (Au, 
Ag, Pt) and dielectric (ZnS) beads can be coated in a controlled way by a silica shell 
[231 12H I2IH I2H E7J- and a dielectric (silica, AU2S) bead can be coated by gold or some 
other noble metal [221I2E1I211E01EI1E2]- One can subsequently etched away silica core of a 
silica-core metal-shell bead and obtain a hollow metallic nanoshell. Either hollow metallic 
nanoshell or a dielectric-core metal-shell bead can be in turn coated in a controlled way by 
the second concentric silica shell (a dielectric overcoat of either metallic shell or a metallic 
core prevents aggregation of the particles by reducing the Van der Waals forces between 
them) [221 EH] and by a further metal shell, thereby forming a "nano-matryoshka" |22j . 
Compared to a simple homogeneous sphere, such a complex multilayered spherical parti- 
cles allows one a lot more freedom in engineering of both elastic [221 12H1 1201 GUB EH] and 
inelastic light-scattering properties [HI]. Many other examples of stratified spheres can 
also be found in nature. As an example, water insoluble aerosols in atmosphere have a thin 
liquid layer adsorbed on their surface. With the addition of an appropriate surfactant, 
water droplets in different hydrophobic solvents (such as oil) form a stable nanometer 
sized structures (with the size depending on the water to surfactant molar ratio), often 
referred to as "reversed micelles" , consisting of a spherical water core coated with a sur- 
factant monolayer [3|. In the case of a biological cell, the appropriate model consists 
of concentric three-layered sphere, corresponding to nucleus, cytoplasm, and membrane 
[To] . The case of a sphere having two coatings is also important for modeling hydrological 
particles coated with biological material and micro-encapsulated material. 

A unique feature of "nano-matryoshka" structures of Prodan et al [22] , which partly 
explains our focus on these structures, is the existence of two coupled nanocavities sur- 
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rounded by metal boundaries (see Fig. [TJ). The latter feature holds promise of large and 
controlled tunability of light-matter interactions, including both the "nano-matryoshka" 
scattering properties and the spectroscopic properties of the atom interacting with such 
a nanostructure. As it will be shown below, such a complex nanoshell geometry strongly 
affects spontaneous emission decay rates, W rad , photostability, and the ratio W rad /W tot , 
known as the fluorescence yield, or simply, quantum efficiency (here W tot = W rad + W nrad 
denotes the total decay rate) (2D]- The spontaneous emission rate W rad is one of three 
Einstein's coefficients (usually denoted as A coefficient). The remaining two Einstein's 
coefficients (usually denoted as B coefficients) describe the stimulated emission probability 
and the absorption probability. When multiplied by the radiant energy density U^duj with 
circular frequency between uo and u + duo, the respective B coefficients then determine the 
stimulated emission and absorption rates at the transition frequency u. Detailed balance 
in thermal equilibrium implies that the knowledge of a single Einstein coefficient is suffi- 
cient to determine the remaining two. Therefore, once the spontaneous emission rate is 
known, the stimulated emission and absorption rates are also unambiguously determined. 
For instance, the stimulated emission probability and the absorption probability are equal 
and the ratio of spontaneous to stimulated emission decay rate remains equal to the mean 
thermal photon number at the transition frequency, a constant which does not change 
with changing environment (see third reference in pQ). Thus, an inhibited (enhanced) 
spontaneous emission necessarily implies inhibited (enhanced) stimulated emission by the 
same factor. There is hope that using complex nanostructures one would be able to tailor 
spontaneous and stimulated decay rates according to one's need and a desired applica- 
tion, such as chemical speciation, LIDAR, fluorescent near-field microscopy, identification 
of biological particles, and monitoring specific cell functions. In this respect, fluorescent 
properties of the atom both inside and outside of a complex nanostructure are of fun- 
damental interest. For instance, by placing fluorescent organic groups or rare earth ions 
(with nm control over the radial position [26:) inside the dielectric core or dielectric shell 
of such a complex nanoparticle (221; a fluorescent nanoprobe can be formed for biophysical 
and biomedical applications [4 j. Alternatively, for the applications in near- field optical 
microscopy |25 : , wherein the probe tip is modeled as a sphere of small radius, one is inter- 
ested in changes in the atom fluorescence properties induced by the presence of a nearby 
complex nanoparticle. 



2 Theory 

In order to characterize the change in the spectroscopic properties of the atom interacting 
with a (complex) spherical scatterer, the frequency shifts, radiative and nonradiative decay 
rates will all be normalized with respect to the radiative decay rates of the same dipole 
but now in a free-space filled in with a homogeneous medium which is identical to that at 
the atom position. Such a normalization of spectroscopic properties brings an advantage 
that, in the case of radiative decay rates, any local field corrections El UH EH] cancel 
out. In the case of a homogeneous dielectric medium (characterized by the refractive 
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index n and dielectric permittivity e) Nienhuis and Alkemade [H7j have showed that 

W r h ad = —W r v a \ (1) 



£ 



where W[ ad and W™ d are radiative rates of the electric-dipole transitions in the dielectric 
medium and in the vacuum, respectively. For a non-magnetic medium, the above relation 
reduces to W rad = nW^ ad . It is emphasized here that these results holds irrespective if a 
homogeneous medium is dispersive or not. (For the general case of a linear, nonconducting, 
absorptive, and dispersive medium see Sec. VIII of Ref. |3~Sj.) 

Formally, irrespective of either the classical or quantum-mechanical descriptions, the 
line broadening and frequency shift of an electric-dipole emitter interacting (via the vac- 
uum electromagnetic field) with a material body can be understood as a result of the 
coupling of the emitted field with its own reflected field. Let us label the concentric shells 
of a complex nanoshell from the nanoshell core outward, with the nanoshell core count- 
ing as the shell number one and the ambient counting as the shell number five. Let rj, 
j = 1, 2, 3, 4, and Sj, fij, kj = UQ^/ejpTj '/ c, j = 1, 2, 3, 4, 5, denote the respective shell radii, 
dielectric permittivities, and wave vectors. Occasionally, as in Fig. the subscript 5 will 
be replaced by h to indicate the host medium (ambient). Then, within the linear response 
formalism of Agarwal and of Wylie and Sipe jHj, the effective shift in the frequency 
separation ujq of two levels of the atom within the n-th shell is given as follows (see Refs. 

0000113), 

ivp' = Re [p ' G(r "' rd ' ' p] = _Re p ' Es(rd ' Uo) ' (2) 

Here r d is the dipole position, p is the transition dipole moment, k n = uj 0yi /e n fi n /c, and 
G(r, r^, id) denotes the scattering Greens function normalized such that the electric field 
E a (r, a>) of the scattered radiation at r due to a dipole p radiating at frequency to at r d is 
given by 

E s (r,uj) = G(r,r d ,£j) • p. (3) 
For a homogeneous sphere one obtains in the quasistatic limit [T3*t IT4] , 

Id — Co>o\ 3 E\ — £2 1 



W r h ad J h 32 e x + e 2 {k 2 r s - k d r d f 

Id — ld \ 3 E\ — 62 1 



W r h ad J ± 16EX + E2 (k 2 r s - k d r d y 



(4) 



where k d is the radiation wave vector in the medium wherefrom the dipole is radiating. 
The most characteristic feature of the quasi-static approximation is that the frequency 
shift is a monotonic function of the dipole distance from the sphere boundary. Moreover, 
for the limiting cases of an atomic dipole in a close proximity to the sphere and in the 
long-wavelength limit the quasi-static approximation fails to account for the retardation 
effects and the radiative decay rate (linewidth) remains unchanged [T31 IT4]. 
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Whereas the change in the effective shift in the frequency separation of two levels is 
given in terms of the real part of G(r, r, u> ), the total decay rate induced by the presence of 
a (multilayered) sphere is determined by the imaginary part of G(r, r, uq). Indeed, within 
the linear response formalism of Agarwal [5] and of Wylie and Sipe jHj, the normalized 
decay rate for the atom within the n-th shell is given as |H1 El EH ECU EE] 

W l 3e 3e 

W h = 1 + 2^k Im [P ' G(l " d ' r "' ^o) • P] = 1 + Im ^3 P ■ E a (r d , u Q ). (5) 

The basic assumption is, of course, that neither the transition matrix element nor the 
transition frequency are appreciably changed by the presence of the interface. 

In the presence of an absorption, as in our case, the decay rate W l comprises the fol- 
lowing two basic decay channels: 1) the process of real (i.e., not virtual) photon emission 
with the photon escaping to the spatial infinity, i.e., radiative decay, and 2) the process 
of real (i.e., not virtual) photon emission accompanied by the subsequent photon absorp- 
tion by the microsphere, i.e., a nonradiative decay, 13 El 13 EH ESI EE]- I n an ideal 
theoretical situation (i.e., a single fluorescent atom, the respective silica and gold shell 
being without any impurities, and multi-photon relaxation absent) the decay rate W t , 
as calculated according to Eq. (jSJ), would be the total spontaneous decay rate. How- 
ever, the Ohmic loss is only one of many other nonradiative mechanisms, such as, for 
instance, multiphoton relaxation, coupling to defects, direct electron-transfer processes, 
and concentration quenching, which may all contribute to the nonradiative decay rate 
W nrad |>jgj ggj gjj g^J but are not included in Eq. ©. Therefore, in practice, W* 
would be the lower limit to the total spontaneous decay rate W tot . 

Obviously, in the absorbing case, the spontaneous decay rate W t , as calculated ac- 
cording to Eq. (JSJ from the imaginary part of Green's function at coinciding arguments, 
G(r, r, oj), pHE] does not coincide with the radiative decay rate W rad . The latter can be, 
up to a proportionality factor, determined as the classically radiated power of a dipole 
which escapes to spatial infinity, or simply the radiative loss, P md ; which a classical dipole 
experiences when interacting with a (multicoated) sphere [TTH ITU IT2"1 ITS] . P rad is calcu- 
lated from the electromagnetic flux given by the surface integral of the Poynting vector 
through a virtual sphere of radius R extending to infinity fUJ El IE] • According to the 
correspondence principle, the radiative decay rates W rad is then given by 

prad 

W rad = . /gx 

niti 

In the absence of absorption, W l and W rad coincide. In the presence of an absorption, 
the ratio W rad /W t , also known as fluorescence yield, is always smaller than one. The 
relative difference between W rad and W t is especially pronounced in the proximity of 
metal boundaries (see Fig. EUl below) . The quantum theoretical expression for the power 
radiated by the spontaneous emission from an excited state in an electric (a magnetic) 
dipole transition is still obtained from the classical expression for the power radiated by an 
electric (a magnetic) dipole, by replacement of the dipole moment by the corresponding 
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transition matrix element. (An expression for the dipole source intensity detected by a 
point detector has been provided by Dung et al [SJ (see Eqs. (34-36) therein).) 

Let ko be the vacuum wave vector and e" (//') be the imaginary part of the dielectric 
function (magnetic permeability) at the observation point. The Ohmic loss contribution 
to the nonradiative decay rates, p nrad ^ is calculated according to formula 



ynrad 



Q(r)dr, (7) 



where the volume integral extends over all the absorbing regions. Q is given as the steady 
(averaged) inflow of energy per unit time and unit volume from the external sources which 
maintain the field, 

Q = ^"|E| 2 + / /|H| 2 ). (8) 

Here the averaging is performed with respect to time and assuming that the amplitude 
of a monochromatic electromagnetic field is a constant. The formula (jHJ) for Ohmic loss 
density remains also valid in the regions of high absorption near resonance frequencies 
of the permittivity and permeability, even when the so-called Brillouin expression for the 
electromagnetic field energy density, 



u = ±. 

8tt 



(9) 



is no longer valid (see Appendix C of Ref. |12j). For simplicity, we will assume that 
/j," = 0, i.e., the Ohmic losses will be entirely determined by an integral of the squared 
amplitude of the electric intensity (see Ref. for calculational details). 



3 Results 

In this section, detailed results of numerical simulations are shown for "nano-matryoshka" 
structures of Prodan et al [22] , i-e., multilayered spheres with silica core and surrounded by 
three additional concentric spherical shells: an inner gold shell, followed by a silica spacer 
layer, and terminated by an outer gold shell. As in Ref. |22]; "nano-matryoshka" struc- 
tures have been considered with the following dimensions: ri/r 2 /r 3 /r 4: : 80/107/135/157 
nm (A), 77/102/141/145 nm (B), and 396/418/654/693 nm (C). As comparative exam- 
ples, the results are also presented for a homogeneous silica sphere with radius r s = 150 
nm (D), a homogeneous gold sphere with radius r s = 693 nm (E), and a homogeneous 
gold sphere with radius r s = 150 nm (F). The radius of the sphere E was chosen to coin- 
cide with that of the sphere C, whereas the radius of the spheres D and F was selected to 
lie between that of the spheres A and B. The radiating wavelength was taken to be 595 
nm, implying gold refractive index ua u ~ 0.248 +i2. 986 [H]. All the spheres are assumed 
to be suspended in an aqueous solution. The respective refractive indices of silica and 
water are assumed to be nsio 2 = 1-45 and nn 2 o = 1-33. There is nothing particular in 
choosing the wavelength of 595 nm, except of it being an emission wavelength of lissamine 
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molecules (13] • Any other (nonresonant) optical wavelength would lead to qualitatively 
similar conclusions. The choice of water as an ambient has been motivated by a fact that 
(i) this often corresponds to an experimental situation and (ii) aqueous solution matches 
biological conditions. 

In numerical simulations, the angular-momentum cutoff value of l max = 60 was used. 
In the case of the total and radiative decay rates, the cutoff value was sufficient to obtain 
convergence on at least 8 significant digits (see Fig. 8 of Ref. J2])- In the case of 
the Ohmic loss contribution, an immediate metal shell proximity provides a numerical 
challenge. However further away from metal interfaces the convergence of up to at least 
8 significant digits can be attained again (see Fig. 9 of Ref. [T2]). 

3.1 Frequency shifts 

The radiative frequency shifts have been calculated directly from the real part of Green's 
function at coinciding arguments according to Eq. (J2|) (see also Eq. (137) of Ref. [T2]). 
The dependence of the frequency shifts on an atomic dipole position inside and outside the 
"nano-matryoshka" structures A-C and the homogeneous spheres D-F for the respective 
radial and tangential dipole orientations is shown in Figs. |2l El Since in the sphere center 
the difference between the radial and tangential orientation of a dipole disappears, the 
corresponding radial and tangential quantities coincide there. 

In order to appreciate changes in the frequency shifts induced by nanoshells A-C 
geometry, let us first discuss the comparative example D of a homogeneous silica sphere. 
Inside the dielectric sphere D, the frequency shift of a radially oriented dipole steadily 
increases from its value of ~ 0.0117 at the sphere origin till that of ~ 261 at r/r s = 
0.995075 (the last sampled point inside the spheres). On the other hand, the relationship 
between the frequency shift and the position of a tangentially oriented atomic dipole has 
an oscillating character, first decreasing from the value of ~ 0.0117 at the sphere origin 
down to ~ 0.00784 at r/r s = 0.388, and subsequently steadily increasing till ~ 133 at 
r/r s = 0.995075. (Though due to the scale of ordinate axis, the shift appears to be a flat 
featureless horizontal line). Outside the dielectric sphere D, it is the frequency shift of the 
radially oriented dipole which exhibits an oscillating behaviour, increasing from the value 
of ~ —334 at the very first sampled point outside the sphere at r/r s = 1.005025, reaching 
the maximum of ~ 0.0038 at r/r s = 1.751294 and then decreasing down to ~ 0.0021 at the 
last sampled point at r/r s = 2.01. For a tangentially oriented atomic dipole, the frequency 
shift gradually increases from the value of ~ —162 at r/r s = 1.005025 up to ~ —0.00017 
at r/r s = 2.01 and compares well with the quasi-static approximation [see Eq. (j3J)]. The 
latter predicts a monotonic increase of frequency shifts from large negative values to zero 
as r/r s increases. Had the sphere radius was larger compared to the emission wavelength, 
one would observe an oscillating relationship between the frequency shift and the position 
of a radiating atomic dipole for all dipole orientations and both inside and outside the 
sphere. In agreement with the quasi-static result for a homogeneous dielectric sphere 
which is optically denser than surrounding medium (e s > eh), the frequency shift of an 
atomic dipole inside the sphere and in close proximity to its boundary is always toward 
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higher frequencies (blue shift) jT3|. On the other hand, the frequency of an atomic dipole 
outside the sphere and in close proximity to its boundary experiences a shift toward lower 
frequencies (red shift) [T3j . 

It is clear from Figs. EJElthat frequency shifts of an atomic dipole interacting with the 
"nano-matryoshka" structures A-C experience significantly enhanced changes than in the 
comparative example of the homogeneous dielectric sphere D. Already at the sphere center 
they can be more than two orders of magnitude larger: pa —1.941 for A, pa —1.691 for 
B, and pa 0.903 for C compared to ~ 0.0117 for D. The emission frequency of an atomic 
dipole in a silica region inside the "nano-matryoshka" structures A and B is always shifted 
toward lower frequencies (red shift). For "nano-matryoshka" structure C, which radius 
is more than four times larger than that of A and B, the red shift is still observed in a 
proximity of gold shells. However, in a marked contrast to "nano-matryoshka" structures 
A and B, further away from gold shells, a small frequency shift toward higher frequencies 
is observed: for r/r s G [0,0.17), (0.69,0.79) and r/r s G [0,0.23), (0.7,0.85) in the case 
of tangentially and radially oriented dipole, respectively. The magnitude of frequency 
shifts substantially depends on the atom position within a dielectric shell. With the 
atom approaching metal shell boundaries, the frequencies exhibit an accelerated decrease 
toward large negative values. We have seen that even in the case of a purely homogeneous 
dielectric microsphere D with a small refractive index contrast, the frequency shifts are 
capable of reaching very high values near the surface of the microsphere. However, for 
the respective "nano-matryoshka" structures A-B, the shifts at the proximity of metal- 
dielectric interfaces can be more than two orders of magnitude larger. 

Outside and at a very close proximity of the outer sphere boundary, large red frequency 
shifts are observed for the complex nanoshells A-C as well as for the homogeneous spheres 
D-F. For instance, for radially oriented dipole source at the very first sampled point 
outside the spheres at r/r s = 1.005025, this red shift ranges from —7392 (B), through 
-5858 (F), -5117 (A), -334 (D), till -67 (C, E). For a tangential dipole orientation, 
the red shift ranges from -3597 (B), through -2839 (F), -2480 (A), -162 (D), till -30 
(C, E). Note that for the largest "nano-matryoshka" C, and a homogeneous gold sphere 
F of the same radius, the comparable frequency shifts at the proximity of the outer sphere 
boundary are 5-times smaller than for the silica microsphere D. Surprisingly enough, the 
quasistatic result (@J) of Klimov et al f3J that, in the close proximity of a spherical 
shells interface, the radial frequency shift is approximately twice as large as the tangential 
frequency shift appears to apply also for complex nanoshells (see Figs. EJHIJ). 

Further away from sphere boundaries, as the value of r/r s increases, the red shift 
typically changes into blue one and vice versa in dampened oscillations around zero. Ac- 
cording to Fig. a noticeably large blue shift with maximum pa 2.14 at r/r s pa 1.154 
is observed for a radially oriented dipole source outside "nano-matryoshka" B, leading 
to repulsive forces between the atom and the dielectric microsphere [H] . The blue shift 
persists in a large interval for r/r s G (1.12, 1.78). For a tangential dipole source orienta- 
tion, this blue shift, which occurs for r/r s G (1.14, 1.275) and (1.63,2.01) with maximum 
pa 0.413 at r/r s pa 1.164, becomes almost five-times smaller. Blue frequency shifts for 
the atom located outside a sphere are also observed in the remaining cases, but they are 
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almost one order of magnitude smaller. For instance, for radially oriented dipole source, 
the excursion above zero does not exceeds ~ 0.148 for C (at r/r s pa 1.144), pa 0.146 for 
E (at r/r s pa 1.154), w 0.059 for A (at r/r s « 1.592), « 0.033 for F (at r/r s pa 1.741), 
~ 0.004 for D (at r/r s « 1.751). In the case of a tangentially oriented dipole source, the 
excursion above zero does not exceeds ~ 0.113 for E (at r/r s pa 1.254), pa 0.104 for C (at 
r/r s w 1.254), « 0.062 for A (at r/r s « 2.01), « 0.05 for F (at r/r s « 2.01), whereas 
for D the frequency shift remains negative till 2.01. Hence, similarly to the case of the 
atom located either close to a plane or inside a homogeneous dielectric spheres [TH], the 
relationship between the frequency shift and the position of the atomic dipole outside the 
spheres considered here has an oscillating character The dampened oscillatory behaviour 
of the frequency shift for the radially oriented dipole outside the sphere D contradicts the 
conclusion reached by Klimov et al |14j . However, Klimov et al jT3] have only studied 
atoms at a distance from a dielectric sphere not larger than r/r s pa 1.2, which is a too 
short distance to observe any oscillating behaviour. Had they drawn frequency shifts for 
larger values of r/r s , they might have observed the oscillatory behaviour, too. 

3.2 Decay rates, the Ohmic loss contribution to nonradiative 
decay rates, and fluorescence yield 

3.2.1 Total decay rate W t 

The normalized total decay rates W l as calculated directly from the imaginary part of 
Green's function at coinciding arguments according to Eq. (jSJ) (see also Eq. (135) of 
Ref. [13)) are displayed in Figs. EH5J Obviously, one finds the rates in the sphere center 
identical for the radial and tangential atomic dipole orientations. An advantage in dealing 
with the normalized decay rates is that any local-field correction El CHI EH] cancels out 
(see also Sec. 14.31 below) and, in principle, a direct comparison between the normalized 
decay rates and experiment can be performed. 

Similarly as in the preceding subsection, in order to appreciate changes in the decay 
rates induced by nanoshells A-C geometry, we will first discuss the comparative example 
D of a homogeneous silica sphere. Inside the dielectric sphere D, the normalized decay 
rate for a radially oriented atomic dipole steadily decreases from its maximum value of 
pa 0.94237 at the sphere origin down to pa 0.83437 at the last sampled point inside the 
sphere at r/r s = 0.995075. On the other hand, the normalized decay rate for a tangentially 
oriented atomic dipole exhibits a weakly oscillating behaviour: it first increases from 
the value of ~ 0.94237 at the sphere origin and reaches its maximum of ~ 0.95173 at 
r/r s = 0.497562, and then decreases down to ~ 0.90179 at r/r s = 0.995075. Outside 
the dielectric sphere D, the normalized decay rate for a radially oriented atomic dipole 
steadily decreases from its maximum value of ~ 1.27798 at the first sampled point outside 
the sphere at r/r s = 1.005025 down to ~ 0.99896 at the last sampled point at r/r s = 2.01. 
On the other hand, the normalized decay rate for a tangentially oriented atomic dipole 
exhibits a weakly oscillating behaviour: it first increases from the value of ~ 0.9824 at 
r/r s = 1.005025, reaches its maximum of ~ 1.00414 at r/r s = 1.860746, and then it 
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decreases down to ~ 1.00366 at the last sampled point outside the sphere at r/r s = 2.01. 
(Though due to the scale of ordinate axis, the decay rates appear to be a flat featureless 
horizontal line in Figs. 01 |3j). Note in passing that for spheres of larger radius compared 
to the emission wavelength one would observe an oscillating behaviour of decay rates for 
any atomic dipole orientation [TUI lllj. 

A characteristic feature of nanoshells is a huge increase of the decay rates for a dipole 
source in a close proximity to metal boundaries, and especially when dipole is within 
the silica core. In a purely dielectric case, such large decay rates are only observed in 
the proximity of sharp resonances of large spheres |10j . whereas in the present case they 
can be achieved with small nanospheres without any special tuning to their internal res- 
onances. For the "nano-matryoshka" structures A and B, the decay rates monotonically 
increase from the values at their center of ~ 0.8751 and ~ 1.7979, respectively, up to the 
respective values at the last sampled core points of ~ 2773 (5412) and ~ 2445 (4765) for 
the tangential (radial) dipole orientation. The normalized decay rates of an atomic dipole 
at the core region of the largest "nano-matryoshka" C exhibit a qualitatively different 
behaviour which is characterized by pronounced minima ~ 0.1696 at r/r s = 0.218955 
and ~ 0.1129 at r/r s = 0.209005 for the radial and tangential dipole orientations, respec- 
tively. Depending on the nanoshell parameters and the atom position, both inhibited and 
enhanced decay rates are observed, with the decay rates maximum values (in hot spots) 
being between two and three orders of magnitude larger than the decay rates minimum 
values (in cold spots) within the same shell. The positional sensitivity of the decay rates 
appears to be more pronounced in the nanoshell core regions than in the second silica 
shell. 

Outside the complex nanoshell C, an oscillatory dependence of the decay rates on the 
dipole position is clearly visible. The behaviour is closely matched by the case E, i.e. the 
case of homogeneous metal sphere of the same radius. The amplitude of the oscillatory 
dependence is much stronger than that discussed earlier for D. A pronounced oscillatory 
dependence of the decay rates on the dipole position outside the sphere can also be seen for 
F and a tangential dipole orientation. For A and a tangential oriented dipole interacting 
with B only a very weak oscillatory dependence is seen. For B and a tangential dipole 
orientation and F and the radial dipole orientation only monotonic decrease of decay rates 
is observed down to w 0.10373 and ~ 0.10188, respectively, at r/r s = 2.01 

3.2.2 Radiative decay rate W rad 

The normalized radiative decay rates W as calculated according to Eq. (jUJ) are shown 
in Figs. |H1 and [7| The radiative decay rate is more pronounced for a radially oriented 
atomic dipole, in which case an order of magnitude enhancement can be expected for an 
optimal atom position inside the second silica shell of the nanoshells A and B. Outside 
the spheres, the largest (an order of magnitude) enhancement of the radiative decay rate 
is achieved in a proximity of the small homogeneous gold microsphere F, closely followed 
by the nanoshell A. The best location of a tangentially oriented atomic dipole appears 
to be within the nanoshell B silica core. Whereas the radially oriented atomic dipole 
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shows typically an enhanced radiative decay rate outside the spheres, in the case of its 
tangentially orientation the radiative decay rate is generally reduced. A strongly reduced 
radiative decay rate is also observed for an arbitrarily orientated atom inside the silica 
shells of the nano-matryoshka structure C. Further away from the spheres outer surfaces 
the radiative decay rate shows dampened oscillations around one for any atomic dipole 
orientation. The radiative decay rate for a tangential dipole orientation exhibits a complex 
behaviour within the nanoshell C silica core, showing there a local minimum of « 0.0204 
at r /r s = 0.199055. 

3.2.3 Ohmic loss contribution to the decay rate 

Compared to the normalized total and radiative decay rates in Figs. |H El El and 
the normalized Ohmic loss contribution to the decay rates shown in Figs. |H1 El does 
not exhibit any oscillatory behaviour outside the spheres. It is reminded here that the 
Ohmic loss contribution to the decay rates has been calculated according to Eqs. (JTj), 
(jEJ). (Further computational details can be found in Sees. 6 and 8.2 of Ref. [P2].) In 
agreement with expectations, the Ohmic loss contribution steadily decreases down to 
zero with the increased atomic dipole distance from the sphere surface. It is also easily 
understandable that the Ohmic loss contribution to the decay rates rapidly increases when 
the atom approaches metal boundary. However, rather contra-intuitively, the Ohmic loss 
contribution within the nanoshell C silica core decreases from its maximum value of 
w 0.2102 at the center down to ~ 0.044 and « 0.047 for the radial and tangential dipole 
orientations, respectively, at the very last sampled silica core point at r/r s = 0.567214. 
For the radial dipole orientation, the Ohmic loss contribution decreases monotonically, 
whereas for a tangential dipole orientation it exhibits an oscillatory dependence with a 
series of local minima and maxima: first decreasing from the central value of ~ 0.2102 
down to « 0.086 at r/r s = 0.228905, increasing up to « 0.137 at r/r s = 0.398060, and 
eventually decreasing down to ~ 0.047 at r/r s = 0.567214. This effect appears to be real 
and not an artifact of computational inaccuracies or numerical instabilities. Calculations 
in extended precision yielded essentially the same result. Note that the local minima of 
the decay rates and the Ohmic loss contribution to the decay rate for a tangential dipole 
orientation occur at more or less the same position within the nanoshell C silica core. 

3.2.4 Fluorescence yield 

In Fig. the ratio W rad /W t is plotted, which is known as the fluorescence yield, or 
simply quantum efficiency. Our W l is that calculated directly from the imaginary part 
of Green's function at coinciding arguments according to Eq. (jSJ) and has earlier been 
shown in Figs. IHE1 However, our W l only comprises the Ohmic loss contribution, W£ rad , 
to the nonradiative decay rate W nrad . The Ohmic loss is only one of many nonradia- 
tive mechanisms, such as, for instance, multiphoton relaxation, coupling to defects, direct 
electron-transfer processes, and concentration quenching, which all contribute to the non- 
radiative decay rate W nrad 03110111111121113 and which may all occur at real experimental 
situations. Therefore, our W t = W rad + W^ Tad does only provide the lower bound for the 
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total decay rate W tot . Since W rad /W tot < W rad /(W rad + W£ rad ), the ratio W rad /W* for 
an averaged dipole orientation plotted in Fig. ITU1 provides an upper theoretical bound on 
the fluorescence yield. Because of absorbing bodies, W£ rad > and the respective ratios 
W rad /W t are always smaller than one. (For the sake of clarity, the trivial non- absorbing 
case D, in which case the ratio is equal to unity for any atom position, has been omitted.) 
As expected, with increasing the atom distance from the outer sphere surface, where the 
Ohmic loss contribution W£ rad decreases to zero, the ratio W rad /W t rapidly approaches 
unity. For all cases considered, the respective fluorescence yields are already larger than 
0.93 at r/r s = 2. At nanoshell centers, a remarkably large value of the fluorescence yield 
(« 0.694) is observed for nanoshell B. In the case of nanoshell C, the fluorescence yield is 
the smallest («s 0.160). As dipole approaches metallic shells, the respective fluorescence 
yields rapidly drop to very small values. This is well-known as fluorescence quenching 
jlj. In spite of decay rates being strongly enhanced in the proximity of metal shells (see 
Figs. El EJ), the decay of the excited atomic states is not accompanied by the emission 
of a real photon, but instead matter quanta are created due to absorption. The fluores- 
cence quenching at the proximity of metal boundaries then implies pronounced maxima 
of the respective fluorescence yields when the atom is located at the middle of the second 
silica shell of nano-matryoshka structures A-C. Indeed, the shell is surrounded by metal 
shells on its both sides. A noticeable feature is also a rather complex behaviour of the 
fluorescence yield within the core of nanoshell C. In the latter case, we have already seen 
a complex behaviour in the case of decay rates and the Ohmic loss contribution to the 
decay rates. The complex behaviour in all these cases is the result of almost 5x larger 
core radius compared to the core radii of nanoshells A and B. 

3.3 Photob leaching and detected intensity enhancement 

If one assumes that photobleaching of a dye takes place only while the dye is in its 
excited states, a sufficiently large enhancement of the spontaneous emission rates can 
significantly lower the probability of switching into nonfluorescent dark (triplet) states, 
thereby increasing stability against photobleaching |3*2~] 150] . The latter means that a 
fluorescent dye molecule can emit more photons before irreversible chemical reactions 
prevent the molecule from any further emission. Let us assume that photostability is 
inversely proportional to the excited state lifetime, i.e., proportional to the total decay 
rate [20 . Then the photostability of the atom interacting with a nanoshell compared to 
that in the free-space filled in with a homogeneous medium which is identical to that at the 
atom position increases by the factor of W tot /W[ ad . On the other hand, there is only the 
probability of W rad /W tot that a given decay will end up with photon at the spatial infinity. 
Given the respective two situations, let N and N be average numbers of the photons 
emitted by the fluorescence source (and detected in the far-field) till photobleaching. Then 
the ratio N/N can be obtained as (W tot /W T h ad ) x (W rad /W tot ) = W rad /W r h ad , i.e., the 
ratio is equal to the normalized radiative decay rate. The latter has been plotted as a 
function of the source position in Figs. |0] and 

The results presented so far have not assumed any incident radiation. An atomic 
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dipole has already been assumed to be in an excited state and the changes in its decay 
properties have been monitored. It is informative to assess the changes in the detected 
fluorescence intensity due to the presence of a metal nanoshell. By taking into account 
that 

1. the respective stimulated emission and absorption probabilities undergo the same 
changes in the presence of material boundaries as the rate of the spontaneous emis- 
sion; 

2. the stimulated emission and absorption rates are given as the product of the inten- 
sity, or squared field strength, with the respective probabilities; 

3. the spontaneous decay rate W l can be enhanced by more than 3 orders of magnitude 
(see Figs. El El), the radiative decay rate W rad an order of magnitude (see Figs. H3 
|7j), and intensity by 4 or more orders of magnitude; 

one finds that an elementary absorption-emission cycle can be accelerated by 7 or more 
orders of magnitude, resulting in the enhancement by 5 or more orders of magnitude of 
the detected fluorescence intensity. 



4 Discussion 

4.1 Size-dependent corrections and nonlocal effects 

The calculations presented so far have been performed assuming the bulk values of gold 
dielectric constant. If theory presented here is to be applied for a multilayered spherical 
particle with a small metallic core or a thin metallic shell with a radius or thickness S < 20 
nm (such as the outer shell of nano-matryoshka B), two effects have to be additionally 
considered. First, the bulk dielectric function is modified, since the electronic mean free 
path is then shorter than in the bulk [33] . Second, nonlocal effects come into play [313 133 
HHj. The first effect can easily be incorporated by replacing the bulk dielectric function 
esiyj) with its size- dependent modification 

2 2 

e{u) = e B {u) + _ x - — ^— — — ■ (10) 

UJ 1 + %UJT B UJ Z + IUJT 1 

Here u p is the bulk plasmon frequency, tb is the relaxation time in the bulk metal, 
t -1 = Tg 1 + vfS^ 1 is the inverse relaxation time (also called damping coefficient T) 
corrected for the finite size of the particle, and vf is the Fermi velocity. More generally, 

T- 1 =T B 1 +AV F S~ 1 , (11) 

where A is a parameter determined by the geometry. For simple Drude theory and 
isotropic scattering one usually takes A — 1. 

On the other hand, nonlocal effects, i.e., when the Fourier transform of the dielectric 
function depends in addition to u also on k, are associated with the resonant excitation 
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of longitudinal bulk plasmon modes (either propagating ones, with frequency above the 
plasma frequency u p , or evanescent ones, with frequency below the plasma frequency u p ) 
[lol WA H5] . The neglect of the nonlocal responses of the substrate is the main reason 
why a phenomenological treatment will generally break down when the radiating atom is 
very close (within a few nanometers) to the sphere surface [IB]. Indeed, the spectroscopic 
properties of the two-level atom interacting with a complex spherical nanoshell have been 
investigated here within the framework of macroscopic Maxwell's equations. The latter 
lead to unphysical results in the limit r — > r s , i.e., when the atom approaches sphere 
surface. For instance, the quasi-static approximation [see Eq. (jl))] predicts that frequency 
shifts increase as (r — r s )~ 3 . The non locality brings about a natural cutoff pole order: 
the nonlocal sphere does not polarize significantly at angular momenta higher than a 
certain cutoff value l c . The latter is of principal importance since it allows for a fully 
converged treatment of multipolar excitation effects Another important feature of 
the nonlocal dielectric function is that it introduces a natural cutoff angular momentum for 
the excitation of near- field modes jlBlES], and thereby a natural cutoff for the convergence 
of the nonradiative decay rates. 

The nonlocal effects may significantly influence both the radiative and nonradiative 
rates of very small nanospheres and nano-matryoshka structures with very thin shells. For 
instance, in the case of a small homogeneous nanosphere with radius of r s < 5 nm, the 
radial and tangential dipolar decay rates for low emission frequencies (ou < 0.5 u p ) can 
be reduced by up to 2 orders of magnitude with respect to the local results [IB]. On the 
other hand, for emission frequencies u > 0.5u p , the decay rates can be up to 2 orders of 
magnitude larger than in the local case [IB]. The theory of Leung [IB] has been extended 
to the case of complex nanoshells (see Sec. 11.1 of Ref. [12]) and its application will be 
presented elsewhere. 



4.2 Radiative decay rates vs the local density of states 

In the ideal case of nonabsorbing and nondispersive material medium (which, however 
(except for the vacuum), does not exist), the trace of the imaginary part of the Green's 
function at coinciding spatial arguments, which enters Eq. (J5J), is related to the local den- 
sity of states, p(u, r) = — (l/7r)Tr Im G(rd, r^, uq), where Tr denotes trace (G is a tensor 
quantity). This relation has prompted claims that radiative decay rates are proportional 
to the local density of states (LDOS). However, even in the ideal nonabsorbing homoge- 
neous case one finds this claim true only when (i) one performs orientational averaging 
over atomic dipole orientations and (ii) properly includes medium dependent prefactors 
[see Eq. (JSJ)]. Once dispersion comes into play, radiative decay rates cease to be propor- 
tional to the LDOS. Indeed, Nienhuis and Alkemade [BZ] have showed that the LDOS in 
the homogeneous dispersive case is 

PM = - 2 ^W), (12) 

where n(u>) = y/ejl is the frequency dependent refractive index of the medium and 
p^(cu) = u 2 /c 3 7i 2 is the LDOS of photon states in vacuum (see Eqs. (29-31) of Ref. 
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[37]). Yet, as it has been emphasized, the Nienhuis and Alkemade formula (JTJ, which 
does not contain any derivatives of the medium refractive index, remains also to be valid 
in the dispersive case. 

The presence of a derivative factor in the LDOS can also be traced down to similar 
derivative factors in the Brillouin expression for the electromagnetic field energy density 
©. Another argument is as follows. Let G(E) = 1/(E — H) denote a scattering Green's 
function, where E is an energy and H is a Hamiltonian. Then, for an isolated eigen- 
value E n of an energy dependent Hamiltonian (our case; see Refs. [S]) the quantity 
— (l/7r)Tr ImG(E) is no longer equal to the Dirac delta function 5{E — E n ), but merely 
proportional to 5{E — E n ) with a prefactor 1/|1 — dH/dE\. The latter is, in general, dif- 
ferent from unity. Therefore, the density of states (DOS) p{E) can no longer be defined as 
p(E) = — (l/7r)Tr ImG(E). Yet another argument is that the integrated density of states 
is given by the expression N(E) = (l/7r)Tr Im In G(E) = — (l/7r)Tr Im \n(E — H). Since, 
for an energy dependent Hamiltonian, the relation (d/dE) [Tr Im In G(E)] = — Tr Im G(E) 
no longer holds, N(E) ^ J E p(E) dE, provided that the density of states (DOS) p(E) is 
defined as p(E) = -(l/n)Tr ImG(E) (see also Ref. [EE]). 



4.3 Local-field corrections 

An atomic dipole couples to the microscopic vacuum fluctuations. If the (unnormalized) 
decay rate W t [see Eq. (J3J] is obtained from macroscopic Maxwell's equation, the differ- 
ence between microscopic and macroscopic vacuum fluctuations can be accounted for by 
local-field corrections. The local field effects always play an important role and constitute 
a critical test of our understanding of the relation between microscopic and macroscopic 
electromagnetic phenomena [21 B]. As discussed by Schuurmans et al .36], one has to 
distinguish between the substitutional and interstitial character of impurities. In the for- 
mer case the well-known empty-cavity factor 3e/(2e + 1) applies, whereas, in the latter 
case, the Lorentz local- field factor (e + 2)/3 is obtained. Since the decay rate W l can be 
expressed in terms of an expectation value of the product of two electric field operators 
(see, e.g., Eq. (25) of Ref. [£]), the local-field factors appear squared. In the near-infrared 
and for optical wavelengths, an inclusion of the corresponding local-field factor is only 
necessary for electric dipole transitions. Indeed, in the above wavelength range, magnetic 
permeabilities p of different materials are all equal to the vacuum value and, therefore, the 
local-field factors for magnetic dipole transitions are trivially equal to one. (For local-field 
corrections in absorbing case see Ref. [TS].) 

However, since we have only discussed normalized decay rate W l [see Eq. (JSJ)], any 
local-field correction cancels out and a direct comparison between the normalized decay 
rates and experiment can, in principle, be performed (see, however, the following subsec- 
tions). 
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4.4 Nonradiative decay rates 

Earlier, at the end of Sec. Ir>.2| it has been discussed that the sum W l = W rad + W£ rad , 
where W™ ad stands for the nonradiative decay rate due to the Ohmic losses, does only pro- 
vide the lower bound on the total decay rate W tot . Consequently, Fig. El only shows the 
upper theoretical limits on the fluorescence yields. Indeed, the Ohmic loss is only one of 
many other nonradiative mechanisms, such as, for instance, multiphoton relaxation, cou- 
pling to defects, direct electron-transfer processes, and concentration quenching, which all 
contribute to the nonradiative decay rate W nrad jSHl EH HH El EH] ■ It turns out that even 
in a purely dielectric case, in the absence of any Ohmic losses and for small fluorescence 
atom concentrations, the nonradiative decay W nrad can be higher than the radiative decay 
W rad , resulting in the fluorescence yield smaller than 0.5 [3J 1221 EH] • When fluorescence 
dye concentrations increase above a certain threshold value, the fluorescence yields of most 
organic dyes are substantially reduced even further with respect to a zero-concentration 
limit value. It should be emphasized that, for a dipole outside the sphere, the radiative 
decay rate W rad is proportional to the intensity of the time-resolved fluorescence spec- 
tra at t — jlH]- Therefore, in the latter case, W rad and W nrad can be disentangled 
experimentally [43J. 

In order to test theoretical predictions experimentally, the choice of a suitable flu- 
orescence source turns out to be a critical issue. Rare-earth ions, which exhibit long 
luminescence lifetimes, are very suitable candidates j2J E2J E3] • However, they are usually 
implanted by an ion deposition resulting in a poor control over their radial distribution 
within a spherical shell. On the other hand, fluorescent organic groups can be placed in- 
side the dielectric core or shell of a complex nanoparticle with nanometer control over the 
radial position j^Ej, but their decay rates are usually strongly affected also by other then 
purely electromagnetic mechanisms, such as concentration quenching. For instance, in 
the case of fluorescein (FITC), the threshold value for a concentration quenching is ~ 0.1 
mM (~1 niM) in liquid (solid) solutions [201 • However pyrene-doped PMMA spheres with 
pyrene concentrations up to 10 mM do not exhibit any concentration dependence |55j . 
Note also that metals per se have inherent photoluminescence j^Zj- However, since metal 
photoluminescence results in a broad band continuum [351 EZ|, it can easily be filtered 
out. 

4.5 Nonlocal decay rates 

The nonradiative mechanisms different from the Ohmic losses do only depend on the 
immediate neighborhood of a radiating dipole. Therefore, one can include all such mech- 
anisms of nonradiative decay rates under the local nonradiative decay rate, W loc . On 
the other hand, the nonradiative decay rate due to Ohmic losses and the radiative decay 
rate can be viewed as a nonlocal decay rate, W nloc = W l = W rad + W£ rad . The rea- 
son behind this nomenclature is that the latter two rates depend on the geometry and 
material composition of the entire sphere and surrounding medium, and not only on the 
immediate proximity of the radiating dipole. The total decay rate is then written as 
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W tot = W loc + W nloc . As a rule, both the nonradiative decay rate due to Ohmic losses 
and the radiative decay rate change if the optical properties of a shell being far away from 
the radiating dipole (for instance, ambient medium) change |J21 E21 EH] • In the case of 
a homogeneous dielectric sphere, the local and nonlocal decay rates can be, in principle, 
disentangled by measuring the total decay rate W tot using the same sphere in different 
environments [21 El E2] ■ For instance, the sphere can be embedded in a refractive index 
matched liquid [21 El E2]- Its radiative decay rates then becomes that of a dipole in a 
homogeneous dielectric slab [531 I56|. Other possible environments include air or liquids 
with different refractive indices [52]. The difference of the respective total decay rates 
measured in different sphere environments then corresponds to the difference of the non- 
local decay rates [22] • The latter is obviously also true for a multicoated sphere. The local 
and nonlocal decay rates can then be separated by fit of Eq. W tot = W loc + W nloc to the 
measured data (121 E2] • Only after the local decay rate W loc is determined, a comparison 
of measured data and theory presented here can be performed [32] . 

4.6 Numerical subtleties 

When dealing with dispersive and absorbing shells, it may happen that the linearly in- 
dependent spherical Bessel functions ji and yi (see Sec. 10 of Ref. [HHj for notations) 
are in fact related by yi(kr) ~ iji(kr) up to almost all significant digits in double preci- 
sion. Consequently, if the spherical Hankel functions hf 1 is later on formed as the sum 
hf\kr) = ji(kr) + iyi(kr) |58j . its precision may be drastically compromised. Therefore, 
it is always recommended to determine h\ f\kr) by a direct independent recurrence, such 
as that proposed by Mackowski et al (SHI [see recurrences (63), (64) therein]. Otherwise the 
radiative decay result for the interior of the nanoshell C may differ by up to four- orders of 
magnitude from the correct one. If one can perform calculations in an extended precision, 
this pathological behaviour can be largely overcome, yet the independent recurrence by 
Mackowski et al [SH] is still highly recommended. 

4.7 Outlook 

Interesting avenue of further research is the study of the effect of sphere's resonances 
on the atomic spectroscopic properties. So far, the latter problem has been thoroughly 
investigated only in the case of a homogeneous dielectric sphere JT3] and a dispersive and 
absorbing sphere (8J. However, a comparative study for complex nanoshells is missing. 
Yet another interesting set of problems arises in connection with strong atom-sphere 
interactions, leading, for instance, to Rabi splitting of an atom levels. Also here only the 
case of a dielectric homogeneous microsphere has been considered |6Uj . 

Another important aspect is the inclusion of nonlocal effects into the treatment of 
spectroscopic properties of the atom interacting with complex nanoshells. As it has been 
discussed earlier, the neglect of the nonlocal responses of the substrate is the main reason 
why a phenomenological treatment will generally break down when the radiating atom 
is very close (within a few nanometers) to the sphere surface [48j . So far, the effect 
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has only been investigated in the case of small metal nanospheres |48j . In this regard, 
complex " nano-matryoshka" structures offer the possibility in studying the nonlocal effects 
for spheres of relatively large radius, provided that metal shells are thin enough. General 
theory for the case of " nonlocal" shells has already been developed [12] and its application 
will be presented elsewhere. 

5 Summary and conclusions 

Frequency shifts, decay rates, the Ohmic loss contribution to the nonradiative decay 
rates, fluorescence yields, and photobleaching of an atomic dipole interacting with "nano- 
matryoshka" structures of Prodan et al have been investigated. The changes in the 
spectroscopic properties of an atomic dipole interacting with complex nanoshells have been 
shown to be significantly enhanced, often more than two orders of magnitude, compared 
to the same atom interacting with a homogeneous dielectric sphere. The detected fluores- 
cence intensity can be then enhanced by 5 or more orders of magnitude. The decay rate 
enhancements can be achieved with small nanospheres without any special tuning to their 
internal resonances. The changes strongly depend on the nanoshell parameters and the 
atom position. Rather contra-intuitively the Ohmic loss contribution to the nonradiative 
decay rates for an atomic dipole within the silica core of larger nanoshells may be decreas- 
ing when the silica core - inner gold shell interface is approached. Surprisingly enough, 
the quasistatic result of Klimov et al fSJ that the radial frequency shift in the close 
proximity of a spherical shell interface is approximately twice as large as the tangential 
frequency shift appears to also apply for complex nanoshells (see Figs. EJ E]). Although 
decay rates are strongly enhanced in the proximity of metal shells, the majority of the 
emitted radiation is absorbed and fluorescence yield exhibits there a well-known quenching 
(see Fig. ITU]) . Although simulation have so far been performed at the wavelength of 595 
nm (the emission wavelength of lissamine molecules 43J), Fig. ITT1 demonstrates that sig- 
nificantly modified spectroscopic properties can be observed in a broad band comprising 
all (nonresonant) optical and near-infrared wavelengths. 

It has been shown here that complex nanoshell structures provide the possibility of a 
controlled tunability in engineering of the radiative decay rates. The ability of controlled 
modification of radiative rates for atoms or molecules in the excited state is of great im- 
portance since dissipative pathways of the excited state can be controlled. For instance, 
one can design nanoprobes with enhanced quantum yield for fluorescent microscopy and 
with enhanced photostability for for biophysical and biomedical applications |U E2] , iden- 
tification of biological particles in fluorescence-activated flow of cytomeres |61j . to monitor 
specific cell functions, or in the cell identification and sorting systems fo^\ Ib3| . Enhanced 
spontaneous emission rates could also provide increased sensitivity in low level fluores- 
cence applications |4U] ISP] . Designing of small noble metal nanoparticles with reduced 
quantum yield, W rad /W tot , at the particle close proximity (cca 1 nm) may have crucial 
implications for the particles use as acceptors in biophysical Forster resonant energy trans- 
fer experiments in vitro as well as in vivo [43J. The theory presented earlier in Ref. ^2] 
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with numerical results shown here may also stimulate design of coated tip geometries for 
applications in near-field optical microscopy [35]. 

The emphasis was on the spontaneous decay, rates. However, detailed balance in 
thermal equilibrium implies that the knowledge of a single Einstein coefficient is sufficient 
to determine the remaining two. Therefore, qualitatively similar behaviour as that shown 
in Figs. El is also expected for the stimulated emission and absorption rates. 

Hopefully, the results presented in this article, in conjunction with computer program 
freely available at http://www.wave-scattering.com, will provide a larger freedom in en- 
gineering of (complex) spherical particles properties, rendering them more suitable for a 
variety of applications. 
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Figure captions 



Figure 1 - A typical spherical complex nanoshell, or "nano-matryoshka" , and its pa- 
rameters. In the present case, the "nano-matryoshka" will be embedded in water and the 
respective shaded and unshaded "nano-matryoshka" concentric regions will represent gold 
and silica shells. 

Figure 2 - The normalized level shifts for the radially oriented atomic dipole radiating 
at wavelength of 595 nm as a function of its distance from the sphere center. The shifts 
are normalized to the radiative decay rate W™ in free-space filled in with the medium at 
the dipole position. For the homogeneous dielectric sphere D, the frequency shift reaches 
the value as large as 261 at r/r s = 0.995075. However, for a better view of the region 
around zero frequency shift, the ordinate axis has been terminated at the frequency shift 
of 10. 

Figure 3 - The same as in Fig. |2]but for tangentially oriented atomic dipole. For the 
homogeneous dielectric sphere D, the frequency shift reaches the value as large as 133 
at r/r s = 0.995075. However, as in Fig. |2l for a better view of the region around zero 
frequency shift, the ordinate axis has been terminated at the frequency shift of 10. 
Figure 4 - Normalized decay rates W t /W[ ad for the radially oriented atomic dipole 
radiating at wavelength of 595 nm as a function of its distance from the sphere center. 
The rates are normalized to the radiative decay rate W™ d of the same atomic dipole in 
the free-space filled in with the medium at the dipole position. 
Figure 5 - The same as in Fig. HJbut for tangentially oriented atomic dipole. 
Figure 6 - The normalized radiative decay rate of the radially oriented fluorescence dipole 
source at wavelength of 595 nm The rate has been normalized to that in free-space filled 
in with the medium at the source position. 

Figure 7 - The same as in Fig. El but for tangentially oriented atomic dipole. 
Figure 8 - The normalized Ohmic loss contribution W^ rad /W[ ad to the nonradiative 
decay rates for radially oriented atomic dipole radiating at wavelength of 595 nm as a 
function of its distance from the sphere center. The Ohmic loss contribution has been 
normalized to the free-space radiative decay rate W™ in the medium identical to that at 
the dipole position. 

Figure 9 - The same as in Fig. |S]but for tangentially oriented atomic dipole. 

Figure 10 - Fluorescence yield (quantum efficiency) at wavelength of 595 nm for an 

averaged dipole orientation as a function of the dipole position. 

Figure 11 - Normalized decay rates W t /W[ ad for an averaged dipole orientation in the 
case of nanoshell A as a function of the dipole position for different emission wavelengths. 
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